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Abstract 

We consider in the article the influence of the irradiation and the 
internal sinks of the point defects on the rate of the flattening of 
the surface structure in solids. The irradiation produces only the 
additional external sources of point defects (vacancies and interstitial 
atoms). The general system of equations is formulated. The solution 
of the system on the stationary stage of the process is found. It is 
shown that depending on the values of some parameters of the solid the 
irradiation can increase or decrease the rate of the surface flattening. 



The reason of smoothing of the surface structure in sohds is the difference 
of the equihbrium point defect concentration on the convex and concave 
parts of the surface. As a result of this difference there are the diffusion 
flows of point defects which tend to redistribute the point defects to attain 
the thermodinamically equilibrium state corresponding the fiat surface. The 
expression for the rate of surface flattening without irradiation was received 
in [1]. Using this expression it is possible to find the surface and volume 
diffusion coefficients analyzing experimentally the rate of surface smoothing 
as a function of parameters of the roughness of the surface [2]. In article 
[3] the irradiation of the solid was introduced into the process of surface 
flattening. The irradiation was considered only as a external source of point 
defects. It was shown in [3] that because the depth of allocation of external 
source of interstitial atoms is bigger then the depth of allocation of external 
source of vacancies the irradiation increase the rate of the surface flattening. 
This conclusion can be changed if we tak;e into account the internal sinlcs of 
point defects. Our article is devoted to this problem: the evolution of the 
surface structure of solids with external sources and internal sinlcs of point 
defects. 

The main system of equations which describes the evolution of surface 
structure was thoroughly analyzed in [3,4]. For small enough intensity of 
irradiation we can disregard the annihilation of point defects (vacancies and 
interstitial atoms) . In this case we can flnd the contribution of vacancies and 
interstitial atoms into the rate of surface flattening separately. Let us flnd 
the contribution of vacancies. 

Following the article [3] we assume that the surface of solid is described 
by the equation: 

^ = /(x,i) = X:^„(i)e— (1) 

n 

where prime in (1) means that n 7^ in the sum (i.e. the equilibrium surface 
is 2; = 0), a; = 27r/A, A is the distance between two identical roughness on the 
surface, 2;-axis is directed into the bulk of the solid. We will consider Znl\ 
as a small parameter and will disregard the second order terms in -2„/A in all 
equations below. Under this assumption the curvature of the surface (1) is 
equal to: 

K{x, t) ^ = ^ Znit)nWe-'''''' (2) 
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To find the diffusion (surface and volume) flows of vacancies which modify 
the surface structure in sohds we should flnd the volume c^ac and the surface 
Uvac vacancy concentrations as the functions of space and time co-ordinate. 
We will find the concentration as the power series in small parameter Zn/ A to 
within the first-order terms in Zn/X. With the surface structure in the form 
(1) such series are 



cU^, z, t) = (z, t)+Y: c^liz, i)e— (4) 

n 

where c!^J^, u!^J^ have zero order in Zn/X and c^^, have the first order in 

Zn/X. 

In was shown in [3] that under real conditions the main system of equa- 
tions describing the spatial and temporal distribution of vacancies Cvac{x, z, t) , 
Uvac{x, t) takes the form: 



dujJc — _,,(',,(0) {^\ _ „, \ < Dv,vac dcUc{z,t) 

z=j{x,t] 

%^ = Dv,vacAc^ac + I{x, Z, t) - 6,ac{Cvac ~ Cq) + C^acV U {x , Z, t) 

Dv,.ac = = fix, t),t)- ^4°) (t)) 

z=f(x,t) 



(Cvac(x, y, Z, t) - ) 



df{x,t) 
dt 



z=f(x,t) 



Ds,vacaAsU^ac + a^^V s{u^acV sU {X, Z, t)) + Dv,vac % 

(5) 

where index "vac" indicates that the quantities correspond to vacancies; 
Dv,vac is the volume diffusion coefficient of vacancies; Dg.vac is the surface 
diffusion coefficient; v is the frequency of vacancy absorption by the surface; 
/ is the intensity of a source of vacancies - the number of vacancies gen- 
erated by irradiation per lattice site and per unit time; 5yac is the rate of 
absorption of vacancies by internal sinks (l/5^ac is the vacancy mean free 
time); cq is the vacancy concentration far from the surface (the equilibrium 
vacancy concentration for the fiat surface); cr = Cq + Co7^ is the equilib- 
rium volume vacancy concentration corresponding to the surface with the 
curvature K; 7 is the surface tension; Q, — a^, a is the lattice constant; k is 



z=f{x,y,t) 
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the Boltzman constant; T is the temperature of the sohd; A5 is the surface 
part of Laplacian A = + dy + d"^; n is the inward normal to the surface 
of the solid; PDv,vac — D' determines the last jump for a vacancies emerging 
from the bulk to the free surface, (5 is the dimensionless coefficient that takes 
into account the presence of a potential barrier for a vacancy emerging at 



df(x,y,t) 



is the contribution 

vac 



the sample surface (0 < /? < 1); = 1/z/; 
of vacancies into the rate of the surface smoothing; U {x, z, t) is the elastic 
potential energy of the vacancies in the elastic field produced by the curved 
surface. Without irradiation the terms containing U {x, z, t) are very small 
but under some conditions for solids under irradiation they can be important. 
The potential energy U {x, z, t) is given by the equation [5]: 

U {x, z, t) = \ ft^ac\pix, z, t) (6) 

where flyac < is the decreasing of the total volume of solid under the 
appearing of one vacancy, p{x, z, t) is the hydrodynamic pressure in solid. 
For isotropic case the pressure p{x, z, t) is the solution of Laplace equation 
[5]: 

Ap{x,z,t)^0 (7) 

with the boundary condition: 

p{x,z = f{x,t))=-fK{x,t) (8) 
Then for the surface of type (1) the potential energy is: 

U{x, z, t) = \n,ac\ 7 E Zn{t)u;''n''e-'^'^e-^^ (9) 

We have used in the system (5) that under the condition 1/ » 
it is possible to take '^nj^Qu!^^{t)exp{—inu;x) — uk — Uq [3], where uk — 
uq + Mo7^ is the equilibrium surface vacancy concentration corresponding 
to the surface with the curvature K. 

We shall find the quasi-stationary solution of the system (5) [3]. We have 
such solution for a long experimental time t » Tmax — max {\'^/Ds,vac, l"^ / Dv,i 
(where I is the characteristic depth at which the source of vacancies is located, 
A is distance between two identical structure on the surface) and for a long 
characteristic time of variation of roughness of the surface (/// >> Tmax) be- 
cause in this case we can neglect the time derivative dcyad dt in the first two 
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equations of system (5) and obtain the quasi-stationary system to within the 
terms ~ Tmax/i- The time t is the parameter at this quasi-stationary stage. 
In other words, we must first determine the quasi-steady-state fluxes for a 
given shape of the surface rehef and then find a closed equation determining 
the time variation of the surface structure on the basis of these fluxes. 
To solve the system of equations (5) wc introduce the new variable: 

i^z-f{x,t) (10) 

and the new function 

V{x, 1 1) = ^) + E ^^"He, t)e-^"'^%„e(x, z, t) - ck(x, t) (11) 

Then the system of equations (5) on the quasi-stationary stage takes the 
form: 



^KKacV') UO,vac) H a M 

AV - f,,V^ - 2f^V,^ = -CK,.. - ^ + -if^ {V + ck- Co) - ^VV^VC/ 



€=0 



D 



V.vac 



9^ 



«=0 



\/(")(a;,^ = 0) = 



df(x,t) 
dt 



= Ds,vacaAsU^ac + a^^V s{u^ac^ sU {X, e, t)) + Dv,vac % 

(12) 

We take the source of vacancies in the form of 5-function: 

/(x,e) = (13) 

where Iq is the number of vacancies emerging at the surface per unit time 
and per unit area. When the sample is bombarded by ions, the profile of 
generated point defects has a Gaussian form, and (7) is a good approximation 
for a narrow distribution of defects. 

We find the solution of the first four equation of the system (12) with the 
accuracy to within 2;„/A [3]. Then we substitute this solution into the fifth 
equation of the system (12) and find the contribution of vacancies into the 
rate of the surface smoothing: 

(14) 



dZn(t) 

dt 
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where qs,vac = A*mc = ^w^n^ + g2^^^ ^nd the condition ag^//? was 

used. 

For interstitial atoms the system of equation (12) has the form [3]: 



AV - f^^V^ - 2/,\4g = 



D 



?=0 

J->V,in fe-i 



D 



= 0) = 



__^^(^V(^\^^0,t)-^u^{t)) (15) 



df{x,t) 
dt 



= a^Ws{Uin^sU{x, C, t)) + Dv,in % 



dV 



^=0 



where we have taken into account that co,m ~ (i.e. co,m << co,t,ac); the 
index "in" in system (15) indicates that the quantities correspond to intersti- 
tial atoms. Solving the system (15) we find the contribution the interstitial 
atoms into the rate of the surface flattening: 



dZn{t) 



dt 



where qs 



(16) 



Then the total rate of surface smoothing is: 



dz„{t) 
dt 



^s,mc^%^^n(t)^'n4e-'-'?-- - DsZ'-f^-Wzn{t)uj^n^e 



2 

s,vac 



-IoZn{t)^l{nin (e~ 



/ H'vac \^ 



V kT 

acQsyVac 



and after integration we find: 



1 z„{t) 
Zn{0) 



t- 



— J-^S,vaca-j^UJ n t — JJv,vac-j^-^ [U U + qg^^ac) 



-/ofi (e"'- 



%nHs,%n 



(17) 

The first and the second terms in equation (17) describe the evolution of 
the surface without radiation. For qs,vac = = qs^in (no internal sinks) these 
terms were received in [1]. The last three terms in (17) is the contribution of 
irradiation into In -^^44. For 0^ vac = Q = Qs in we rederive the result of article 
[3]. 
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For qs,,ac « un {qs^in « ton), q^^vac « l/4ac (g^.m << l/kn) and 
tolvac << 1 {tolin « 1) we find from tlie equation (17) tlie corrections to tlie 
results of [1] and [3]: 



1 ^n{t) Co-fQ 4 4 CojQ 33/1, llvac \ . 

/qI^ 7 I r^^jac I 4 4/_. , . „ /oi^7l^m| 4 4/-1 T \, 
--L'S,vac 77f, to n [i- — tvacqs,vac)t' — L>S,in 77^'^ ^ [i- — Hnqs,in)t — 

T r\ 2 2 \ 1 /i Q.s,in , Qs,in \ ; f Qs,vac , Qs,vac \ \ , /io\ 
- IqUoJ n Uin\l h - ; ^ - Ivac 1 h _ ' „ > t (18) 



We see that we have the second order correction ql^ac/^"^ ^'^ result 
without sinks and without irradiation (the second term in equation (18) [1]) 
and the first order correction qs,vac/<^, lvac(ls,vac {qs,in/to, hnQsM ) to the result 
with irradiation and without sinks [3]. 

The sign of the last term in eq. (17) depends on the sign of the expression: 

For ip > the last term is negative ant it increases the rate of the surface 
smoothing and for < the last term increases the roughness of the surface. 
The leading term in the third term have a;^-dependence on cu. 
For qs,vac » oj {qs,in >> <^) the equation (17) takes the form: 

1 Zn{t) Co7l^ 4 4 Co-fVl 2 2+ 

VJO) " ~^S,vaca-j^OJ n t - Dv,vac-j^(ls,vacOJ U t- 

- Io^^^^Y~ {hne'^'""^'-'" - /mce~'''"=*^''"'=} t (20) 

where the conditions Ivac^^'^ << Qs,vac and Ivac^^^ « <ls,vac were assumed. In 
this case the second and the last terms in equation (19) have the same uP' 
dependence on the frequency uj. The sign of the last term is determined by 
the expression: 

^1 = /me~'*"^^'*" - lyacd'^"'"'''"'"""' (21) 
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Depending on the values of qs^in {qs,vac) and /«„ (Ivac) the irradiation tends 
to decrease (for ■i/'i > 0) or to increase (for ipi < 0) the roughness of the 
surface. The expression (15) shows that the interstitial atoms tends to make 
the surface become flat and the vacancies created by irradiation tends to 
increase the roughness of the surface. The influence of the point defects on 
the evolution of the surface is determined by the value of /e^''^^.For example, 
if the interstitial atoms are absorbed by internal sinks stronger that the 
vacancies are then the vacancies are the main defects and the expression ipi 
is negative - the irradiation increases the roughness of the surface. 

In conclusion, we derived the expression for the rate of the surface smooth- 
ing taking into account the external sources of point defects (due to irradi- 
ation) and the internal sinks of point defects. We have found that the ir- 
radiation has tow contribution into the rate of surface smoothing: the first 
contribution (the third and the fourth terms in equation (17)) increase the 
rate of the surface flattening and the influence of the second contribution of 
irradiation (the last term in equation (17)) on the evolution of the surface 

structure depends on the sign of the expression -0 = lin ^1 — -|- ^i^^ — 

" ^ + (^O^ ^^'Vac « qs,in « and qs,vaJvac « 1, 

qs,vaclvac << 1) and the expression Vi = /j„e~'*"«^''" - Cce~''""=^'''''"= (for 
Qs.vac » qs,in » ■ FoT i/j > {ipi > 0) thc sccoud contribution of 
irradiation increases the rate of the surface flattening and for -0 < (-01 < 0) 
the it increases the roughness of the surface. 
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